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Abstract 

The temporal dynamics of double ionization of H2 has been investigated both experimentally and 
theoretically with few-cycle laser pulses. The main observables are the proton spectra associated 
to the + H"*" fragmentation channel. The model is based on the time-dependent Schrodinger 
equation and treats on the same level the electronic and nuclear coordinates. Therefore it allows 
to follow the ultrafast nuclear dynamics as a function of the laser pulse duration, carrier-envelope 
phase offset and peak intensity. We mainly report results in the sequential double ionization regime 
above 2 x 10^^ Wcm~^. The proton spectra are shifted to higher energies as the pulse duration 
is reduced from 40 fs down to 10 fs. The good agreement between the model predictions and 
the experimental data at 10 fs permits a theoretical study with pulse durations down to a few 
femtoseconds. We demonstrate the very fast nuclear dynamics of the Hg" ion for a pulse duration 
as short as 1 fs between the two ionization events giving respectively H2" from H2 and + H+ 
from H2'. Carrier-envelope phase offset only plays a significant role for pulse durations shorter 
than 4 fs. At 10 fs, the laser intensity dependence of the proton spectra is fairly well reproduced 
by the model. 

PACS numbers: 33.80.Rv, 33.80.Eh, 42.50.Vk 
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I. INTRODUCTION 



Recent outstanding advances in ultrafast laser physics have led to the generation of few- 
cycle pulses in the near-infrared and visible ranges [l|, l2j and in the XUV range 3|, |j] . Very 
efficient spectral techniques are used for a complete reconstruction of the electric field in 
the temporal domain. In the infrared and visible ranges, these techniques are known as 
frequency-resolved optical gating (FROG) ^] and spectral phase interferometry for direct 
electric field reconstruction (SPIDER) |6|. In the XUV range, the conventional optical 
elements used for FROG and SPIDER cannot be used. In consequence, the attosecond 
pulse is reconstructed from spectral informations based on the photoionization of atoms in 
the gas phase. In the first attosecond experiments, the spectral phase was measured through 
two-photon, two-color photoionization of atoms 0, 4|. More recently, an extension of the 
FROG concept to the XUV photoionization was proposed and successfully applied to the 
measurement of attosecond pulses 7|, |8| . 

Following these developments, this report deals with the temporal dynamics of double 
ionization of the H2 molecule induced by near-infrared few-cycle intense laser pulses in 
the 10^^-10^^ Wcm-2 intensity range. The strong field response of the H2 molecule and 
its related molecular ion and isotopic species have been extensively studied in the past 
jsl, |lO, 11|. However, the H2-laser interaction remains a subject of great interest for few- 
cycle laser pulses because of the attosecond electronic and femtosecond nuclear time scales 
and the relatively simple decay channels (ionization and dissociation) leading to H^, + H, 
and 11+ + H"^. For instance in D2, the D+ + D channel was used to propose the idea of an 
attosecond molecular clock based on the rescattering dynamics of the first ionized electron 
leading to the dissociation of the D^ ion 12|]. The same excitation scheme and dissociation 
channel were recently studied in order to control the electron localization with 5-fs carrier- 



envelope-phase-locked laser pulses 



13j. We have recently shown that the proton spectrum 



allows to detect the presence of a pre-pulse or a post-pulse using a 10-fs pump-probe 
scheme [l^ and the important concept of charge resonance enhanced ionization introduced 
by Bandrauk et al. (l5 |. 

In H2, double ionization leads to two bare protons. In consequence, our main experimental 
and theoretical diagnostic is the proton spectrum as a function of the pulse duration, carrier- 
envelope offset phase, and intensity. Double ionization induced by few-cycle laser pulses 
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exhibit two regimes : nonsequential double ionization at low laser intensities be 



ow 2 X 



10^'* Wcm ^ and sequential ionization at higher laser intensities 

m ll8|, ll9|, 120, |2l|. 

Charge resonance enhanced ionization belongs to the sequential ionization regime. However, 
it does not occur for few-cycle laser pulses because of the ultrashort pulse duration. After 
the first ionization event of H2, the resulting molecular ion does not have enough time 
to stretch and to reach the internuclear distance range where enhanced ionization takes 
place [15 1 . Nonsequential double ionization has been studied in detail with an emphasis 

n □ 

on the rescattering dynamics [17|, 119|. Here we propose an experimental and theoretical 
study of the sequential regime at intensities above 2 x 10^'* Wcm~^ where rescattering is 
less important 18|]. The interesting feature is the time delay between the first and second 
electron removals. For few-cycle laser pulses this time delay might be expected to be a few 
femtoseconds corresponding to the pulse risetime. In the meantime, a nuclear wave packet 
arises in the H^ molecular ion from the nonresonant coupling of the ground electronic state 
X^S+ and the first excited dissociative state A^E^. Therefore the delayed ejection of the 
second electron takes place during the evolution of the molecular ion internuclear distance. 
The resulting proton spectrum is shifted to lower energies in comparison to what might be 
expected from an instantaneous two-electron ejection. 

In spite of the simplicity of the above picture , a quantitative prediction of the proton 
spectrum as a function of the ultrashort pulse parameters is not straightforward because of 
the complicated nonlinear couplings leading to single and double ionizations. In addition, 
the evolution of the nuclear wave packet of the intermediate H^ ion has to be included in the 



theoretical framework. Usually this problem is solved using a two-step approach |16l . Il8 |. 
After the first ionization event, the nuclear evolution of the H^ ion is solved by numerical 
integration of the time-dependent Schrodinger equation. The initial wave packet is given by 
the projection of the H2 ground vibrational state onto the different ion vibrational states. 
The second ionization event leading to double ionization is then calculated as a function 
of the time-dependent vibrational wave packet. Tong and Lin give a clear account of the 
procedure in Ref. ISj. Here we propose a complementary approach for the quantitative 
analysis of our experimental results. Double ionization is treated using a unified approach 
based on a two-electron model where the internuclear distance remains a full quantum 
variable in order to extract ^ the nuclear dynamics during the interaction. This model has 
been used in a recent paper [22] for the study of the mechanisms leading to double ionization 
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with 1-fs laser pulses. This manuscript focuses on other effects and on a more detailed 
comparison with experiments. 

The paper is organized as follows. The experimental set-up is presented in Section [IT] 
including the 10-fs pulses generation set-up and the time-of-fiight detection of the ions. 
The theoretical model is described in Section IIIII with an emphasis on a realistic field-free 
molecular description. In Section IIVI the experimental results are compared to the theoret- 
ical predictions for a better understanding of the ionization and fragmentation mechanisms. 
In Section HV] pulse duration, carrier-envelope phase, and intensity dependences of the pro- 
ton spectra are successively presented and discussed. Since the carrier-envelope phase of 
our 10-fs pulses is not controlled experimentally, this dependence will be commented from 
theoretical predictions only. The conclusions are finally summarized in the last section. 



II. EXPERIMENTAL SET UP 



A. Laser system and pulse compression method 

The ultrashort pulse generation set-up is based on a 1-kHz titanium:sapphire laser system 
and a hollow fibre pulse compression stage. The 1-kHz laser chain is built following the 



conventional chirp-pulse-amplification scheme [23|]. It consists of an oscillator, a stretcher, 
a regenerative amplifier, and a compressor. The system delivers pulses with energies up to 
600 yuJ, duration of 40 fs, and a central wavelength of 795 nm. The pulse compression stage 
is designed following techniques introduced by Nisoli et al. ^ and Sartania et al. The 
laser beam is focused with a 700-mm-focal-length lens onto the tip of a 700-mm-long hollow 
fibre with a 250-/im inner diameter. The hollow fibre is housed on a V-groove in a chamber 
filled with argon gas. The non-linear Kerr effect in argon leads to self-phase modulation and 
wavelength spectrum broadening while the fibre waveguide ensures a spatially-homogeneous 
spectral broadening. Optimum argon operating pressures were found around 700 mbar for 
40-fs and 600-/iJ input laser pulses. 

Figure [1] shows the wavelength spectra before and after the hollow fibre filled with argon. 
Self-phase modulation leadsJ;o a noticeable broadening. The usual multi-peak structure is 

2|] . The oscillatory behavior comes from interferences of waves 



observed as in other studies 







with the same frequency but generated at different times within the laser pulse. The overall 
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FIG. 1: (Color online) Wavelength spectra before (dashed curve) and after the hollow fibre (full 
curve) filled with argon gas. The argon pressure is 700 mbar. The laser pulse energy and duration 
are respectively 600 /iJ and 40 fs. 

transmission of the hollow fibre set-up including all the optical elements is above 50 %. 
After recollimation by an / = 1 m concave silver mirror, pulses are recompressed to 10 fs 
using several reflections on commercial broadband chirped mirrors. The pulse duration is 
measured using a home-made interferometric autocorrelator. Figure [2] presents an interfer- 
ometric autocorrelation of 10-fs pulses corresponding to the wavelength spectrum presented 
in Fig. [TJ The measured autocorrelation signal is compared to a calculated interferometric 
autocorrelation assuming a constant spectral phase. The good agreement between the mea- 
sured and calculated curves shows that the chirped-mirrors compression stage works well for 
the second-order group delay dispersion compensation. The remaining disagreement comes 
from cubic and quartic residual phase that cannot be compensated with our set-up. The 
laser pulse duration relative uncertainty is estimated to be ±20 %. Ultrashort 10-fs pulses 
with energies above 200 fiJ are therefore available for subsequent experiments. 

The pulse duration measurement presented in Fig. [2] has been recorded after the opti- 
mization of the group delay dispersion. This optimization is only valid for the autocorrelator 
because the chirped mirrors compressor is tuned in order to compensate for the air travel 
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FIG. 2: (Color online) Interferometric autocorrelation (lAC) of 10-fs laser pulses. The experimental 
and calculated signals are given respectively by the full and dotted curves. The lAC curve is 
calculated using the Fourier transform of the frequency domain electric field, assuming a constant 
spectral phase. 

to the autocorrelator and also for the small amount of group delay dispersion introduced 
by this device. The proton spectra are recorded in a vacuum chamber which is located 
elsewhere in the laboratory and which introduces a different group delay dispersion than 
the autocorrelator. Therefore the pulse duration will have to be optimized in situ at the 
location where the laser-molecule interaction takes place. The method is to introduce more 
negative group delay dispersion than necessary and then to compensate for it with a variable 
thickness of fused silica which exhibits a positive dispersion of 36.1 fs^rad~^mm^^. Proton 
spectra are systematically recorded for different fused silica thicknesses. As expected, the 
overall proton spectrum is shifted to the highest energy for the shortest pulse duration. This 
optimization procedure will be commented in more detail in Section [IV] since it involves a 
thorough understanding of the molecular response. 
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FIG. 3: (Color online) Covariance map from H2 recorded with linearly-polarized 40~fs laser pulse 
at / = 2 X lO^'' Wcm~^ and p(H2) = 5 x 10^® mbar. The collection electric field is 25 Vcm~^. 
The bottom and left panels represent the usual proton time-of-flight spectrum and the central 
panel represents the covariance map. The fl, f2, and f3 peak labels correspond to protons ejected 
towards the detector. The bl, b2, and b3 labels are associated to protons ejected backwards the 
detector. 



B. Ion detection set— up 



The ultrashort pulses are sent in an ultrahigh vacuum chamber equipped with a 75-mm- 
focal-length on-axis parabolic mirror which allows to get laser intensities up to 10^^ Wcm~^. 
The hydrogen gas is introduced through an effusive gas jet at very low pressures down to 
3 X 10^^'^ mbar which is the residual pressure of the chamber. Molecular hydrogen ions and 



protons are detected using a 1150-mm-long time-of-fiight spectrometer based on the Wi 
McLaren configuration and devoted to experimental studies of multiple ionization |2J, 



25|. 



Fragmentation channels and the associated kinetic energy release spectra are determined 
using covariance mapping introduced by Frasinski et al. 26|, |27|]. In particular for H2, it 
is important to separate the + H~*" double ionization channel from the + H single 
ionization channel. 
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We recall here the main features of this technique. First of all, it allows to work with 
more than one ionization event per laser shot. Fragments coming from the same dissociation 
channel are expected to produce time-of-flight signals that fluctuate following a same pat- 
tern on a shot-to-shot basis. For fragmentation channels involving two detected bodies, the 
covariance coefficient C2{Ti,T2) therefore measures the statistical correlation of fragments 
arriving at times Ti and T2. The method is illustrated in Fig. [31 The conventional time- 
of-flight proton spectrum is represented in the bottom panel. The time-of-flight signal has 
been delayed by 3500 ns in order to record only the signiflcant proton signals. The proton 
spectrum is highly symmetric around the relative time of flight 220 ns and exhibits 3 peaks 
of protons emitted towards the detector and labeled fl, f2, and f3. The peaks labeled bl, 
b2, and b3 are associated to protons emitted backwards the detector. Let us recall here that 
the proton time of flight T is given by T = Tq ± P/eFc, where Tq is the time of flight of a 
proton with a zero initial momentum, P is the modulus of the initial momentum, the ± sign 
is positive for a proton emitted away from the detector and negative for a proton emitted 
towards the detector, Fc is the collection electric fleld {Fc = 25 Vcm~^ in these experiments) 



and e is the elementary charge 27|. In addition our spectrometer exhibits a strong angular 
discrimination due to its large longitudinal dimension. Therefore it allows the detection of 
protons with an initial momentum parallel to its axis. Since the angular discrimination is 
stronger for backwards protons, the backwards peaks bl, b2, and b3 in Fig. [3] are slightly 
smaller than the fl, f2, and f3 peaks associated to forwards protons. 

The covariance coefficient C2{Ti,T2) is represented in the central map in Fig. [3l For 
an easy visualization of the correlations, the conventional proton time-of-flight spectrum is 
also represented in the vertical left panel. A correlation between peaks f3 and b3 clearly 
appears in Fig. [3] and corresponds to the Coulomb explosion channel H"*" + H+. The peaks 
fl and f2 and their associated backwards peaks bl and b2 are not correlated with any other 
peaks. They are associated to the H+ + H dissociation channel. Following the relationship 
between the time of flight and the momentum, the proton kinetic energy release spectra are 
extracted from the time-of-flight data using the associated time-to-energy normalization 
factors. Some proton energy spectra are represented in Fig. HJ The bottom curves represent 
the total and the + H'^ covariance energy spectra extracted from the time-of-flight data 
presented in Fig. [31 Although there remains some covariance noise for protons coming from 
+ H, the proton spectrum associated to the H"^ + H+ channel is unambiguously identifled 
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FIG. 4: (Color online) Proton spectra from H2 recorded with linearly-polarized 40-fs laser pulse 
at / = 2 X 10^^ Wcm~^ and / = 2 x 10^^ Wcm~^. For each laser intensity, the full curve represents 
the total proton spectrum and the dashed curve represents the covariance proton spectrum from 
the H'^ + dissociation channel. 

as the broad peak above 1.5 eV of the conventional proton energy spectrum. Therefore in 
the following, we only present conventional proton spectra which exhibit a better statistics 
than the covariance spectra. 

III. THEORETICAL TREATMENT 

This section introduces the theoretical background and numerical methods we have de- 
veloped for studying strong field Coulomb explosion dynamics of the hydrogen molecule. 
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A. The molecular system and the laser— molecule interaction 



We follow the double ionization dynamics of the H2 molecule by solving the time- 
dependent Schrodinger equation for the electronic and nuclear motions 

ih ^X.*(i^, r, t) = [Ho + Mnt] Xs'^iR, r, t) , (1) 

where Ho is the field-free Hamiltonian and Vint the field-molecule interaction potential. 
The body-fixed coordinate r = {ri, refers here to the electrons, while R represents the 
internuclear vector. The electronuclear wave packet is denoted by Xs^iR^f^t), where Xs is 
an antisymmetric two-electron spin wave function, while the spatial wave function ^(.R, r, t) 
is symmetric with respect to the exchange of the two electrons (singlet state). This approach 
does not assume a separation of the electronic and nuclear coordinates, thus going beyond 
the usual Born-Oppenheimer approximation. 

The Hamiltonian of the field-free diatomic molecule is expressed as 

no = fn + ^ + Vi2 (r-i - r2) + J] , (2) 



i=l,2 



where T„ is the nuclear kinetic operator and V^iri — the inter-electronic repulsion. The 
mono-electronic Hamiltonians hi are expressed as the sum of the electronic kinetic operator 
and the electron-nuclei interaction potential 

hi^-—W?+VeniR,ri). (3) 
zm 

The ionization dynamics is initiated by the length gauge radiative coupling 

^nt{rur2,t)^-e{rl + r2)■E{t) (4) 

associated, in the dipole approximation, with the linearly polarized classical electric field 

E{t) = Eof{t)cos{Lut + (p)e, (5) 

where Eq is the field amplitude, cu the angular frequency, and (p the carrier-envelope offset 
phase. The pulse shape is given by the Gaussian-like expression 

m = sin^ , (6) 
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with total pulse duration r. The frequency u corresponds to a central wavelength of 800 nm 
and the internuclear coordinate R is constrained along the field polarization vector e. The 
two electrons, of coordinates rj = also assumed to oscillate along this axis. Recent 

numerical studies have indeed shown that two-electron dynamics in molecules is mainly 



characterized by a one-dimensional motion in linear polarization [28|, |29 |. 

In order to mimic the dynamics of the real H2 molecule, we have introduced two -R- 
dependent softening parameters a{R) and l3{R) in the Coulomb potentials describing the 
electron-electron 

Vuizi - Z2) = [{zi - Z2f + a\R)\ (7) 
and electron-nuclei interactions 

VUR: Zi) = -Y, [(-2^ + ^ X ^/2)' + (^\R)\ • (8) 

s=±l 

These two softening parameters are assumed to vary slowly with the internuclear distance. 
In a first step, the parameter [3{R) is adjusted at each internuclear distance in order to 
reproduce the energy of the ground electronic state of [30]. The parameter a{R) is 
then obtained by reproducing the energy of the ground electronic state of the hydrogen 
molecule jSlj] . These potential curves, presented in Fig. [5l confirm that the adjustment of 
these two parameters allows for an accurate reproduction of the exact molecular potentials 
even though the electronic problem is presently reduced to a single dimension. The variation 
of these two parameters with the internuclear distance is given in Fig. El 

B. The electronuclear wave packet propagation and the initial H2 wavefunction 

In order to calculate the single and double ionization of the hydrogen molecule submitted 
to an intense and pulsed laser radiation, we propagate the total wave function \E'(-R, r, t) in 



time during the entire pulse using the split operator method developed by Feit et al. [32 1 



where the total Hamiltonian Hq + Vi^t = T + V{t) is split in two parts corresponding to the 
kinetic and potential propagations 

xe-'^''l^'' + o{5t'). (10) 
12 




FIG. 5: (Color online) Calculated potential curves of the electronic ground states of H2 and H^. 
Squares are obtained from the model presented in the paper and full curves are the ab initio 



energies reported in Refs. 



30, 



3l| 



T represents here the sum of the nuclear and electronic kinetic energy operators, while Vit) 
includes all potential operators. The kinetic and potential propagations are performed in 
the momentum and coordinate spaces respectively. Three-dimensional Fast Fourier Trans- 
formation (FFT) allows rapid passage back and forth from one representation to the other 
at each time step. Typical grids extend up to z"^^^ = z^^^ = 100 ao and = 10 ao with 
(512)'^ grid points. A time step of 5t ^ 1 as is necessary to achieve convergence. For very 
short pulse durations, an additional field free propagation is performed after the end of the 
pulse in order to give enough time for the ionized electrons to reach the asymptotic region 
were the wave function can be analyzed (see Sec. lIII Cl hereafter) . 

The initial wave function \l/(i?, 7",^ = 0) is taken as the ground nuclear and electronic 
state of the hydrogen molecule, calculated within the present non-Born-Oppenheimer model 
using the imaginary time relaxation technique [33|] . The equilibrium internuclear distance of 
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FIG. 6: (Color online) Soft-Coulomb parameters. The full curve and the dashed curves represent 
respectively the electron-nuclei and the electron-electron smoothing parameters. 
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FIG. 7: (Color online) Probability density from the vibrational ground state of H2. Full curve: 
This model. Squares: Calculated from the ab initio ground electronic state u = of H2 given in 



Ref. 



3li. 
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H2 , -Re — 1.4 ao, is perfectly reproduced, and Fig. [7] shows the very good agreement obtained 
between the nuclear probability density 



P{R) = JJ \^{R,r,t = 0)f dzidz2, (11) 

calculated here and the probability density xl{R) of the ground vibrational state v = 
of Hg calculated from the Born-Oppenheimer ground electronic state potential given in 
Ref. m. 



C. The wave packet analysis 

The single and double ionization probabilities are analyzed using a well-established car- 
tography technique js^ . The plane {zi,Z2) is partitioned in three regions Fq, Fi and F2 
corresponding respectively to H2, H^, and Hg"*". Double ionization occurs in the asymp- 
totic region F2 = {\zi\, \z2\ > za} with za = 20ao ^ Re- The single ionization region Fi 
is defined as {\zi\ < za, \zj\ > za}, and the neutral H2 molecule is found in the region 
Fo = 1-221 ^ -2a}. a schematic illustration of these three regions is given in Fig. [8]which 
represents the partition of the (^1,^2) plane. The outgoing ionization flux is accumulated 
during the time propagation in the regions Fi and F2 to extract the single and double ion- 
ization probabilities as a function of time. Absorbing boundaries are imposed at the end of 
the grid to avoid spurious reflection effects. 

To extract the proton kinetic energy distributions S{E, t) obtained after Coulomb explo- 
sion, we use a simple mapping which relates S{E, t) to the probability density 

P2{R,t) = jj^ \^{R,r,t)\^ dzidz2, (12) 

in the F2 region, using the Coulomb relation E = 0.5/ R and the requirement of particle 
conservation 

P2{R,t)dR = S{E,t)dE. (13) 

The kinetic energy E denotes here the energy of a single proton. The energy distribution 
S{E,t) is finally accumulated over the entire time propagation to obtain the kinetic energy 
release spectrum S{E) which is measured in experiments. 

Figure [9] represents the comparison between an experimental proton spectrum recorded 
at 4 X 10^^ Wcm~^ with a pulse duration of 10 fs and two spectra calculated at the same laser 
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FIG. 8: (Color online) Partition zones of the electrons coordinates z\ and z^- The three regions 
To, Ti and r2 correspond respectively to the H2, + e~, and + e~ + e~ systems. 

intensity with pulse durations of respectively 10 fs and 12 fs. The proton peak associated 
to the + channel is better reproduced with a 12-fs pulse duration than with a 10-fs 
duration. It is quite noticeable that a relatively small variation of 2 fs leads to a measurable 
proton energy shift. Concerning the agreement with the experimental data, the very high 
sensitivity of the proton spectra might explain the discrepancy between the experimental 
and theoretical proton spectra at 10 fs when one takes into account the relative accuracy of 
the pulse duration and laser intensity measurements. In the following, we choose to compare 
the experimental data with calculations performed with a pulse duration of 10 fs in order to 
deal with a theoretical model without any adjustable parameter. The agreement at different 
laser intensities will be commented more thoroughly in Section IIVI 



IV. RESULTS AND DISCUSSION 



A. Pulse duration dependence of proton spectra 

This section is devoted to the dependence of the proton spectra as a function of the pulse 
duration. The comparisons between different pulse durations are performed while keeping 



16 



the same peak laser intensity. For instance, Fig. [TO] represents two proton spectra recorded 
at 8 X 10^^ Wcm~^ with respectively 40-fs and 10-fs laser pulses. Both spectra exhibit the 
above-mentioned separation between the + H and H"*" + H"*" dissociation channels. In 
particular for this last channel, the proton spectra have maxima at 3.8 eV at 40 fs and 5.5 eV 
at 10 fs. We first would like to emphasize the large shift of 1.7 eV when the pulse duration 
is reduced from 40 fs to 10 fs. It is also important to notice that double ionization cannot 
be considered as an instantaneous double ionization process at 10 fs since such a sudden 
electrons removal would produce a proton spectrum peaked at 9.2 eV. This instantaneous 
two-electron emission would lead to the proton spectrum represented in Fig. [TT] by the far 
right curve. The corresponding calculation is based on the projection of the population of 
the f = vibrational state of H2 represented in Fig.[7]onto the H+ + repulsion curve with 
the appropriate normalization factors. The interpretation is here straightforward: After the 
first ionization step, enough time is left to the ion for a significant stretching before the 
second ionization event. The proton peak shifts towards higher energies at 10 fs since the 
internuclear distance is reduced because of the shorter pulse duration. 

To the best of our knowledge, intense infrared laser pulses with durations below a few 
femtoseconds are not yet available. Theoretical predictions thus become highly desirable 




Proton kinetic energy (eV) 

FIG. 9: (Color online) Proton spectra from H2 at 4 x 10^^ Wcm~^. Full curve: Experimental with 
a pulse duration of 10 fs. Dotted curve: Theoretical with a pulse duration of 10 fs. Dashed curve: 
Theoretical with a pulse duration of 12 fs. 



17 



in order to know the H2 molecule sensitivity to ultrashort pulses. Figure [TT] shows five 
calculated proton spectra from right to left which correspond to the instantaneous two- 
electron ionization for the far right spectrum and then to pulse durations of respectively 
1 fs, 2 fs, 4 fs, and 10 fs. The laser peak intensity remains 8 x 10^^ Wcm"^ as in the 
experimental results in Fig. [TUl All curves have been normalized to unity for an easier 
comparison. The calculations are performed with a zero carrier-envelope offset phase. The 
most striking feature comes from the proton spectrum calculated with a 1-fs pulse duration. 
The spectrum peaks at 8 eV and is therefore already shifted by 1.2 eV towards lower energies 
in comparison with the proton spectrum from instantaneous double ionization. Even for such 
an ultrashort laser pulse, nuclear motion takes place and leads to a measurable shift of the 
proton spectrum. As for the experimental spectra presented in Fig. [TUl the theoretical 
proton spectra are shifted to lower energies as the pulse duration is increased from 1 fs to 
10 fs. The spectrum calculated et 10 fs peaks around 5.3 eV, again in good agreement with 
the measurement shown in Fig. [TOl 
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FIG. 10: (Color online) Proton spectra from H2 at 8 x 10^'' Wcm ^ recorded with pulse durations 
of 40 fs and 10 fs. 
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FIG. 11: (Color online) Proton spectra from H2 calculated at 8 x 10^"^ Wcm~^ for different pulse 
durations. From right to left: "0", 1, 2, 4, and 10 fs. The "0"-fs pulse duration spectrum represents 
what is expected from an instantaneous two-electron ionization of H2. 

The above experimental and theoretical results give an a posteriori justification of our 
experimental method of the pulse optimization. Indeed, the optimization of the ultrashort 
pulse duration can be performed in situ inside the ion spectrometer using the proton spec- 
trum since it shifts to higher energies as the pulse duration is decreased. After the initial 
pulse duration optimization using the interferometric auto correlator, the proton spectra are 
recorded with different thicknesses of fused silica inserted in the beam just before the ion 
spectrometer entrance window. An optimum of the proton spectrum to higher energies is 
systematically looked for. If such an optimum is not found, then additional negative group 
delay dispersion is introduced in the compressor by adding additional bounces onto the 
chirped mirrors. Figure [12] represents the results of the optimization procedure using fused 
silica plates with thicknesses up to 5 mm. The maximum shift of the proton spectra is 
found with a thickness of 3 mm and corresponds to the shortest pulse duration available 
within our set-up inside the ion spectrometer. Moreover, one can infer from Fig. [12] that 
the total number of detected protons presents a minimum at the shortest pulse duration 
for a fused silica thickness of 3 mm. Indeed the shortest pulse duration is associated to a 
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FIG. 12: (Color online) Optimization of the ultrashort pulse duration. Proton spectra from H2 
recorded at 2 x 10^^ Wcm~-^ with different thicknesses of fused silica from to 5 mm inserted in 
the beam path before the ion spectrometer. 

reduced internuclear range and hence to a higher energy gap in order to reach the double 
ionization H"^ + H"^ threshold. In addition, charge resonance enhanced ionization plays 
a minor role in comparison with larger internuclear distances produced with longer pulse 
durations. Finally, such a procedure was used with other molecules such as N2 and O2 and 
represents a straightforward method for pulse duration optimization [25|]. 



B. Carrier envelope phase dependence 



The carrier-envelope phase dependence of strong field effects was first investigated in 
high-order harmonic generation and above-threshold ionization 35|, l36|] . Concerning molec- 
ular dissociation, Roudnev et al. have theoretically shown that it is possible to control the 
and HD+ dissociation with the carrier-envelope phase 37|. More recently a first exper- 
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FIG. 13: (Color online) Calculated proton spectra from H2 at 8 X 10^^ Wcm-2 function of the 
carrier-envelope phase i-p and the pulse duration. Full curves: 99 = 0. Dashed curves: i-p = tx 12. The 
left (a), middle (b), and right (c) top panels correspond to different pulse durations of respectively 
1 fs, 2 fs, and 4 fs. The vertical arrows mark the positions of the different peaks. The corresponding 
electric field is represented in the lower panel (d) as a function of time for a pulse duration of 1 fs. 

imental evidence was given in the dissociative ionization of D2 [13]. The carrier-envelope 
phase dependence of the intramolecular electronic motion leads to the localization of the 
electron on one or the other nucleus. Therefore, the resulting dissociation channels + D 
and D + D+, where the left or right position of the deuteron D"*" indicates its initial emis- 
sion direction, can now be separated as a function of the carrier-envelope phase. Finally, 
Tong and Lin investigated the carrier-envelope phase dependence of nonsequential double 



ionization of H2 by few-cycle laser pulses 38|. They found that the strong dependence of 



the double ionization yields is due to the return energy of the rescattering electron. 

Here we address the question of the carrier-envelope dependence of the proton spectra 
from double ionization at laser intensities where double ionization is mainly sequential. In 
the nonsequential regime, Tong and Lin have shown that the proton spectra lie within the 



same proton energy range [38[]. Moreover at laser intensities above 1.5 x 10^^ Wcm~^ and 
a pulse duration of 5 fs, proton spectra from nonsequential double ionization are found 
to be independent on the initial phase. We theoretically confirm this tendency at higher 
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FIG. 14: (Color online) Proton spectra from H2 at different laser intensities 8 x 10^^, 4 x 10^^, 
2 X 10^^, and 10^^ Wcm~^ from top to bottom. Full curves: Experimental. Dashed curves: 
Theoretical. 

laser intensities in the sequential regime. Figure [T3] represents proton spectra calculated at 
8 X 10^^ Wcm~^ for two values of the carrier-envelope phase (p = and (f = 7i/2 rad. In 
addition, calculations were performed for three values of the laser pulse duration (a) 1 fs, 
(b) 2 fs, and (c) 4 fs. A significant shift of 1.2 eV is only observed for the 1-fs pulse whereas 
no energy shift is visible for a pulse duration of 4 fs. In addition a smaller production of 
H"*" ions is calculated with = 7r/2 as compared with = for 1-fs and 2-fs pulses. One 
can also notice in Fig. [TST d) which represents the time dependence of the associated electric 
fields in the case of the 1-fs pulse, that a higher peak intensity is achieved with ip = 0. 
In consequence, the sequential ionization of the ion is more delayed at (f = tt/2 than 
at (f = 0, thus resulting in a higher proton energy when ip = 0. However as the number 
of cycles increases within the pulse duration, this effect disappears and the proton energy 
range is no more dependent on the carrier-envelope phase. 
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C. Laser intensity dependence 



Double ionization of hydrogen by infrared laser light is a highly nonlinear process. The 
intensity dependence therefore follows a complicated behaviour which has to be studied 
both from combined experimental and theoretical works. Figure [T3] presents four experi- 
mental and theoretical proton spectra from 10^^ to 8 x 10^^ Wcm~^. The spectra are shifted 
towards higher energies as the intensity is increased. The experimental spectra are max- 
imum at 4.2 eV, 4.9 eV, 5.0 eV, and 5.5 eV at respectively 10^^, 2 x 10^^, 4 x 10^^ and 
8 X 10^^ Wcm~^. This effect was also observed in multiple ionization of nitrogen and oxygen 



with pulse durations of 40 fs and 10 fs and comes from sequential electron emission [25|, |39 |. 
The time delay between the first ionization and the second ionization, leading respectively 
to and to H+ + H+, is reduced because the necessary instantaneous intensity for the 
second ionization comes sooner when the laser intensity is increased. 

The experimental spectra are reasonably well reproduced by the model predictions with 
some differences mainly in the high energy sides of the spectra. In spite of these discrep- 
ancies, the model predicts the correct energy range of the expected proton spectra at 10 fs 
for different laser intensities. It therefore constitutes a worthy predictive tool for the use of 
ultrashort pulses when one considers the very complicated dynamics of the hydrogen molec- 
ular response. Several effects might explain the observed differences in addition to the laser 
and ion measurements uncertainties. The calculation is performed for one laser intensity 
and thus does not take into account the intensity distribution within the focal volume. This 
might explain the discrepancies in the high-energy wings of the spectra where nonsequential 
double ionization does play a noticeable role at low laser intensity 2^. Indeed for a given 
peak laser intensity, low intensities are distributed over a much larger volume than high 
intensities and may noticeably contribute to the observed discrepancies in Fig. [Ml 



V. SUMMARY AND CONCLUDING REMARKS 



We have presented a detailed experimental and theoretical analysis of H2 sequential double 
ionization using ultrafast laser pulses in the 10^^-10^^ Wcm~^ intensity range. We have 
shown that the resulting proton spectra are very sensitive to the temporal, phase, and 
intensity characteristics of the pulse. More precisely, higher energy protons are emitted 
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when the pulse duration decreases and when the peak intensity increases. On the other 
hand, the effect of the carrier-envelope phase offset is only significant for pulse durations 
shorter than 4 fs. 

The main characteristics of the measured proton spectra from double ionization of H2 
are well reproduced by a full quantum calculation based on the time-dependent Schrodinger 
equation. The i?-dependent soft-Coulomb parameters introduced in this study make the 
model relatively simple and manageable using tractable computing facilities. These param- 
eters are fitted using ab initio calculations of the ground electronic states of H2 and as a 
function of the internuclear distance R. We are confident that the good agreement between 
experiment and theory at 10 fs allows to give valuable predictions for shorter laser pulses 
and in particular about the ultrafast dynamics of the ion as a function of the pulse 
duration, carrier-envelope phase offset and peak intensity. 

The degree of control of few-cycle laser pulses becomes more and more sophisticated 
[isl ]. In this respect, the H2 response to the associated high laser field might become a 
complementary diagnostic to ultrashort laser pulses in addition to the existing FROG and 
SPIDER techniques. In particular the relative simplicity of the proton detection associated 
to a reliable model might be of some interest in ultrafast laser physics. 
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